Negative differential resistance in normal narrow bands - superconducting junctions 

with Andreev reflection 
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We have calculated current-voltage characteristic curves for normal-superconducting junctions 
with Andreev reflections and different types of electronic bands. We found that when the normal 
band is narrow, of the order of the superconducting energy gap, a negative differential resistance 
appears at a voltage of the order of the band width plus the gap value. In case two bands contribute 
to the total current the conductance can be smaller than unity at voltages above the gap value. Our 
simulations may provide an answer to different experimental data of the literature that were not 
yet understood. 
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The measurement of the current-voltage (J — U) char- 
acteristic curves of junctions made from similar or dif- 
ferent materials separated by an insulating layer is one 
of the usual transport methods to obtain information 
on the electronic band structures near the Fermi level 
Ep. If one of the materials of the junctions is a su- 
perconductor, the I — U curve provides a direct mea- 
surement of the energy gap A, for example. In this 
case the voltage region smaller than A/e (e is the elec- 
tronic charge) the conductance G is in general much 
smaller than Gn, the value obtained at U 3> A/e. In 
case one has a plain normal-superconducting junction 
without any insulating intermediate layer, the conduc- 
tance G(U < A/e) - 2G(U > A/e) = 2G N , a rather 
surprising phenomenon taking into account that for a 
normal-normal metals junctions the normalized conduc- 
tance G/Gn = 1 at all voltages. 

Andreev studied theoretically the heat transport 
through a normal-superconducting plain junction and ob- 
tained the factors 2 and 1 for the normalized conductance 
in these two voltage regions Blonder, Tinkham and 
Klapwijk (BTK) [2j studied the reflection of an electron 
at this kind of junctions and obtained a similar solution 
for the conductance, i.e. it decreases from a factor 2 to 
1 when the voltage grows from zero to values above the 
energy gap of the superconducting part of the junction. 
The physical explanation for the value of 2 at U < A/e 
is related to the reflection and transmission of the elec- 
tron of the normal part into the superconducting one. 
The electron current from the normal part transforms in 
a Cooper pair current in the superconductor and a hole 
reflects back into the normal part. Using a different ap- 
proach Garcia, Flores and Guinea Q studied a similar 
phenomenon observed in scanning tunneling microscopy 
and obtained the same equation as in Ref. [2j. 

All the I — U curves discussed above were obtained 
for the case of ballistic transmittance between a normal 
metal and a superconductor. This approach has been ex- 



tended for the case of ferromagnets with a certain spin 
polarization and for a diffusive regime in a recent publica- 
tion [4. In the last years there have been several reports 
[EL WvA lil. lil. [ToL 11, 12 1 showing experimental data with 
conductance values less than its normal state value Gn 
at U > A/e. Although some modeling for this behav- 
ior on the basis of non- ideal interfaces, proximity effects 
and energy losses were proposed, the generality of this 
behavior observed in different kind of samples suggests 
that other explanation may be needed. 

Assuming two different main cases for the normal- 
superconducting junctions, see Fig. [TJ we calculate the 
I — U characteristic curves and conductance as a function 
of bias voltage U. Our aim is to show in which cases the 
conductance gets smaller than unity (in reduced scale, 
i.e. G/Gn) at voltage values above the gap implying the 
contribution of a negative differential resistance to the 
total current. We show below that this can be observed 
when the band width W of the normal material is of the 
order of the energy gap of the superconducting material, 
what we call a normal narrow band. The conductance 
exhibits a negative differential resistance with a maxi- 
mum absolute value at A < ell < A + W, specially in 
the two dimensional case and it tends to zero for voltages 
eU > W + A. 

We assume that the junction is composed by two differ- 
ent materials one of them a superconductor with a energy 
gap A. Following Refs. 
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the Andreev conductance at 
zero temperature (in units of quantum of conductance) 
G(U) is given by 



G(U) = 1 - 
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which depends on the transmittivity T in the normal 
state; a = (U/A) - [{U/A) 2 - l] - 5 . 

The current J as a function of the voltage U at the 
junction can be calculated taking into account that it is 
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FIG. 1: Energy band structures at the normal (N) - super- 
conducting (S) junctions. We used two main models for the 
normal part of the junction: (a) Usual N-S junction, (b) Junc- 
tion between a superconductor and a semiconducting material 
with a band width 2W. Both graphs are valid at T — K 
and at no applied bias voltage U. 
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controlled by Cooper pairs Andreev currents if eU < A 
However at higher U the current is controlled by 
quasiparticles and for U > 7 A / e we have practically con- 
duction between normal materials 0, 0, y[ • The solution 
for the current I is 



FIG. 2: (a) Theoretical conductance G vs. normalized volt- 
age eU /A for T — 1 calculated using Eq. ([TJ (the case of 
Fig. [Ha)), (b) The same as in (a) but for T = 0.5 and nor- 
malized to its value at eU S> A. (c) Total current I calculated 
with © for T = 1 and (d) T = 0.5. 
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where A is a constant that depends on the junction ge- 
ometry and on the integration average on a, the angle 
that the incoming particles form to the interface and 
B = N(0)evF, with N(Q) the density of states at Ep and 
vf the Fermi velocity. This integral corresponds to bal- 
listic 3D case for a system in which the occupied normal 
band is much broader than A, i.e. Iw^A- The transmit- 
tivity can be expressed as T = 1/(1 + Zq), where Zq was 
defined in Q. Notice that Zq = and oo means T = 1 
and 0. 

For the case of Fig. QJa) the conductance (calculated 
assuming a transmittivity T = 1 , see Eq. (JTJ) has the 
known value of 2 at U < A and reduces gradually to 1 
at voltages above A, see Fig. [U[a). The total current / 
through the junction is calculated following Eq. and 
is depicted in Fig. HJc). In case we have T — 0.5 we get 
the results shown in Fig. [2jb,d) for the conductance and 
current. With exception of the curves in Fig. 6, all others 
are normalized to their normal-state values. 

In the case we have a normal metal with a narrow 
band width W ~ A, then the situation changes and the 
equation for the current is given by 



Am* 



(W-eU') b (E F +eU') b G(U')dU' , (3) 



where the parameter b — 0.5, 0, —0.5 corresponds to the 
3D, 2D and ID case and the parameter A has the same 
meaning as in Eq. ^ and m* is the effective mass of 
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FIG. 3: Continuous lines: Total current / vs. normalized 
voltage eU /A for the case of a normal part with a narrow 
band W = 2A for T = 1 and 0.5 for the 3D (a,b) and 2D 
case (c,d). Dashed lines represent the calculated conductance 
dl/dU. These calculations correspond to the case of Fig.QJb). 



the carriers. The correction term (Ep + eU') b can be in 
general neglected if Ep ^> eU. The theoretical results 
below are calculated at zero temperature. 

If we have a narrow band in energy with its minima 
at U = near the superconducting gap, the integral to 
estimate the current is like in Eq. but with the bottom 
of the band that depends on the assumed band width 
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FIG. 4: (a) Total current I vs. normalized voltage eU/A 
for T = 1 for the case of a normal part having two materi- 
als touching in parallel the superconducting part, one with 
W ^> A and a second with W = 1.5 A, calculated using 
I = Iw»a + 0.5Jw/~a (Eq. J2} and Eq. Q) for the 3D case, 
(b) The same as in (a) but for T — 0.5. (c) and (d) show the 
corresponding conductances dl/dU. In (d) the conductance 
is normalized by its value at ell 3> A. 



FIG. 5: (a) Total current I vs. normalized voltage ell/ A for 
the case of a normal part having two materials touching in 
parallel the superconducting part, one with W 3> A and a 
second one with W = 1.5 A for T = exp(— 1) for the 2D case, 
(b) The same as in (a) but for T — exp(— 0.5). (c) and (d) 
show the corresponding conductances dl/dU normalized by 
its value at elJ ^> A. 



2W and therefore of the order of A. Note that in this 
case the current will not be symmetric for positive and 
negative voltages. For positive voltages Fig. [3] shows the 
current for the case of a conductor with a narrow band 
(see Fig. W(b)) filled up to W = 2A in the 3D and 2D 
cases for two different transmittivities. Note that in both 
3D and 2D cases there appears a negative differential 
resistance region at 2 < eU /A < 3 for T = 1 and at 
slightly higher voltage values for T = 0.5, see Fig. [3] 
The behavior of the conductance dl/dU depends on the 
transmittivity value as well as the dimensionality of the 
normal band, see dashed lines in Fig. [3J In case the 
voltage is negative we will obtain the same curve if the 
Fermi energy were at the middle of the band. Otherwise 
we need to do the calculations for the specific case. 

One could also treat the case in which the narrow band 
material is present in parallel to the larger band one, 
both contributing in parallel to the current through the 
junction. This is an interesting case that may occur at 
the interfaces of the contacts. When the two materials, 
the normal and superconducting one, are put together 
narrow bands at the Fermi energy may appear upon ma- 
terials used, junction geometry and quality. The current 
and the conductance are presented in Figures 0] to [6] for 
different cases. 

In the particular cases treated below we have reduced 
the transmittivity of the narrow band respect to the large 
one by a factor of two. This assumption means that the 



effective mass of the carriers in the smaller band is half of 
that in the larger one. Nevertheless, we will see that the 
negative differential resistance contribution from the nar- 
row band influences the total conductance. We note that 
qualitatively similar characteristic curves as we describe 
below were observed in many experiments but in general 
their meaning was not discussed. The cases we discuss 
below can really happen because these narrow bands can 
be formed at the contacts between the normal and su- 
perconducting parts. Figure 2] shows the calculated total 
current / given by the sum of two normal conductors in 
parallel, one with a large band width W A and the 
other with W = 1.5A, for T = 1 (a) and T = 0.5 (b). 
The figures (c) and (d) show the calculated conductances 
dl/dU. Figure [5] shows the results for the 2D case and 
for T = exp(-l) or T = exp(-0.5). 

We describe now the last case when the small and 
large bands contribute to the current but their carri- 
ers have to overcome a barrier of height <j> existing be- 
tween them and the superconductor. If <fr is not very 
large then it will allow more current as thepotential U 
increases. An example can be seen in Ref. [12j. Figure [6] 
shows the results of the calculations for two bands using 
T = exp(-(2 - OAU) - 5 ). It can be seen that at low 
voltages the conductance is smaller than at high volt- 
ages with an intermediate region where it decreases due 
to the negative resistance contribution from the narrow 
band. This case is different from the previous cases, see 
for example Figs. IHand^l because there the bands just 
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FIG. 6: Total current vs. normalized voltage for a normal 
conducting part with two parallel contributions with W 3> A 
and W = 1.5A with a transmissivity T = exp(-(2-0.4*7) a5 ) 
for the 3D case (a) and 2D case (b). The dashed line in (a) is 
the conductance G following Eq. {l}. The figures in (c) and 
(d) show the corresponding conductances dl/dll. 



dicates that our model should be useful to understand 
experimental results. 

For wide band ferromagnetic materials and ballis- 
tic transport (which produces minima for bias voltages 
smaller than the energy gap) further extensions of the 
BTK model were reported in Ref. approaches that 

can be further developed including the contributions of 
narrow bands. A narrow band effect can be also treated 
taking into account a proximity effect [10(. However, if 
the conductance minima are at higher voltages than the 
equivalent to the energy gap (as e.g. the case of Fe/Ta 
[5j), the proximity effect can unlikely explain it, but it 
is a clear sign for a narrow band effect. These possibil- 
ities have not been yet discussed in any of the existing 
experimental papers and stress the possible existence of 
interesting physics at contacts and their interfaces. 
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contribute with different transmittivities but these do not 
depend on U, i.e. there is no tunnel barrier that change 
the conductivity appreciably. The increase of the con- 
ductance at large U values, as is the case of Fig. [51 may 
occur because when a contact is formed it could have 
some oxide that acts as a potential barrier. 

In conclusion, in this paper we have discussed the 
contribution of different types of bands of normal- 
superconducting junctions with Andreev reflections. The 
results show: (1) For a single normal wide band W 3> A 
the conductance G and the current vs. applied voltage 
behave as expected, i.e. G = 2 for ell < A and tends to 
1 for eU > A. (2) For a junction with a normal narrow- 
band width 2W ~ A and at eU > W + A the current 
tends to a constant and the conductance to zero. At in- 
termediate voltages we find a region of negative differen- 
tial resistance (current decreases with voltage) as shown 
in Fig. [3] (3) As two bands, one with a large and the 
other with a narrow band width, contribute in parallel, 
a variety of cases appear. In any of those discussed here 
the negative differential resistance part coming from the 
narrow band has a clear influence on the expected I — U 
characteristics. (4) Finally, if we assume a not very high 
potential barrier between the normal and superconduct- 
ing parts, the conductance does not saturate but steadily 
increases at high voltages eU > W + A. These results 
indicate that two electronic bands at the Fermi level can 
complicate the form of the characteristic I — U curves 
in real junctions. A comparison of our results with the 
available experimental data 0, 0, 0, H, H, [13, U , IH m - 
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